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Abstract. We study competition between the dissipative and coherent effects in 
the entanglement dynamics of two qubits. The coherent interactions are needed for 
designing logic gate operations with systems like ion traps, semicondutor quantum 
dots and atoms. We show that the interactions lead to a phenomenon of periodic 
disentanglement and entanglement between the qubits. The disentanglement is 
primarily caused by environmental perturbations. The qubits are seen to remain 
disentangled for a finite time before getting entangled again. We find that the 
phenomenon is generic and occurs for wide variety of models of the environment. We 
present analytical results for the time dependence of concurrence for all the models. 
The periodic disentanglement and entanglement behavior is seen to be precursor to 
the sudden death of entanglement (ESD) and can happen, for environments which do 
not show ESD for noninteracting qubits. Further we also find that this phenomenon 
can even lead to delayed death of entanglement for correlated environments. 
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1. Introduction 

It is now well understood tliat entanglement is the key resource for implementation of 
many quantum information protocols like teleportation, cryptography, logic operations 
and quantum communications [H El El HJ [5l [6] . Bi-partite entanglement i.e entangle- 
ment among two quantum mechanical systems each envisaged as a quantum bit (a 
quantum mechanical two level system analogous to a classical bit), has been found to 
be particularly important in this context. Numerous methods of producing qubit-qubit 
entanglement have been investigated during the past decade. A method, which is of 
particular interest in the context of quantum logic gate operations with systems like 
ion-traps and semiconductor nanostructures, relies on the coherent interactions among 
the qubits [9l[l0l[7l[8l[IIl[T2l[T3l[Tl[l5l[l6]. An earlier proposal by Barenco et. al. [7] 
has shown how one can implement a fundamental quantum gate like the C-NOT gate 
using dipole-dipole interaction among two quantum dots modeled as two qubits. This 
was followed by another proposal from DiVincenzo and Loss [8J in which they showed 
how the Heisenberg exchange interaction between two quantum dots can be used to 
implement universal one and two-qubit quantum gates. In their model the qubit is real- 
ized as the spin of the excess electron on a single-electron quantum dot. They proposed 
the electrical gating of the tunneling barrier between neighbouring quantum dots to 
creat an Heisenberg coupling between the dots. Finally they showed explicitly how by 
controlling the exchange coupling one can implement a quantum swap gate and XOR 
operation. Moreover they also showed the implementation of single qubit rotation using 
pulsed magnetic field. Further in a later work Cirac and ZoUer [9] discovered that by us- 
ing the coulombic interaction among two ions one can implement a two-qubit quantum 
logic gate operation. Clearly many proposals require interacting qubits for two qubit 
quantum gates. 

However for a computation to progress efficiently one needs sustained entanglement 
among the qubits as they dynamically evolve in time. This can be achieved effectively 
if the quantum mechanical system under evolution is weakly interacting with its sur- 
rounding. In practice though as the system evolves the system - environment interaction 
becomes stronger thereby inhibiting loses in its initial coherence. This loss of quantum 
coherence is known as decoherence [17] and leads to degradation of entanglement. Thus 
the study of dynamical evolution of two entangled qubits coupled to environmental de- 
grees of freedom is of fundamental importance in quantum information sciences. In 
recent years numerous studies have been done in this respect [HI [191 1201 [211 [221 [23l [2i] . 
One study in particular predicted a remarkable new behavior in the entanglement dy- 
namics of a bi-partite system. It reported that a mixed state of an initially entangled two 
qubit system, under the influence of a pure dissipative environment becomes completely 
disentangled in a finite time [23]. This was termed as Entanglement Sudden Death 
(ESD) [25] and was recently observed in two elegantly designed experiments with pho- 
tonic qubits ^26j and atomic ensemble [27j. Note that an earlier proposal have discussed 
pausible experiment to observe ESD in cavity QED and trapped ion systems [2H|- The 
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phenomenon of ESD have motivated numerous theoretical investigation in other bipar- 
tite systems involving pairs of atomic, photonic, and spin qubits [291 1301 [31], [32] , multiple 
qubits [33] and spin chains [HI [35l [36] • Further ESD has also been studied for different 
environments including collective vacuum noise [38], classical noise [39] and thermal 
noise [IHl [21 [12] . Moreover random matrix environments have been studied [131 [B] . 
These authors [H], also point out the differences in the time evolution of concurrence 
arising from the internal dynamics of two entangled qubits due to the level splitting of 
each qubit. 

ESD in continuous variable systems has also been extensively studied. In particular 
the problem of oscillators interacting with different environments has attracted lots of 
interest [22l [13 [13 [13 [HI [H] . Note that the conditions leading to ESD and probable 
ways of suppressing it are currently being actively investigated [28l [50]. In particu- 
lar it has been shown how ESD can be avoided by using external modulation with an 
electromagnetic field [511 [521 [53] and can even lead to sudden birth for some cases [51] . 
Moreover sudden birth of entanglement has also been predicted for structured heat baths 
[551 [56] and certain choice of initial conditions of the entangled qubits [57]. In another 
recent work it has been shown that under a pure dephasing environment for a general 
two mode N-photon state ESD does not occur [58]. This result was explicitly proven 
for a general 3-photon state of the form |^) = a|30) + 6|21) + c|12) + d\03). 

Even though numerous investigations on ESD in a variety of systems have been 
done so far, the question of ESD in interacting qubits remains yet open. In this paper 
we investigate this question for a system of interacting qubits in contact with various 
models of the environment. We show that due to coherent qubit-qubit interactions two 
initially entangled qubits, get repeatedly disentangled and entangled as they dynami- 
cally evolve leading to dark and bright periods in entanglement [59]. Moreover we find 
that the amplitude of bright periods reduce with time and eventually at some finite 
time vanishes completely, thereby causing death of entanglement. Our investigations 
also reveal that the length of the dark periods depends on the initial condition of the 
entangled qubits and also on the interaction strength. Further we find dark and bright 
periods in entanglement in presence of interaction among the qubits, for initial states 
which do not exhibit sudden death but simple asymptotic decay of entanglement in 
absence of the interaction. We find the existence of dark and bright periods to be generic 
for interacting qubits and occurs for a wide variety of models for the environment. We 
show this explicitly by considering various models of the environment which induce cor- 
related decays, pure and correlated dephasing of the qubits. All of these models exhibit 
the phenomenon of dark and bright periods even though some of them don't show ESD. 

The organization of this paper is as follows. In section II we discuss the model 
for two interacting qubits in contact with a simple dissipative environment and formu- 
late their dynamical evolution by solving the quantum-Louiville equation of motion. In 
section III we develope the theory to study the dynamics of entanglement of the two in- 
teracting qubits and calculate the time evolution of the concurrence under the influence 
of environmental perturbations. In section IV we then study the entanglement dynam- 
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ics of two interacting qubits under the influence of pure dephasing environment. We 
find that coherent qubit-qubit interaction not only leads to dark and bright periods in 
entanglement it also delays the onset of ESD. Further in section V we do a detail study 
of the dynamics of qubit-qubit entanglement for both non-interacting and interacting 
qubits for two different correlated models of the environment. In section V A we focus 
on dissipative environments inducing correlated decay of the qubits. Here we find that 
for non-interacting qubits there is no ESD and even though entanglement vanishes for 
certain initial conditions at some instant, it gets partially regenerated quickly and then 
decays very slowly. When we include the interaction among the qubits we find that 
entanglement exhibits the phenomenon of dark and bright periods. We further study 
the behavior of two qubit entanglement for a pure correlated dephasing environment in 
section V B. We find that the correlated dephasing leads to delay of ESD in absence of 
qubit-qubit interactions. We see that the degree of delay depends on the strength of the 
correlation. Here again when we include the qubit-qubit interaction we observe dark 
and bright periods in entanglement with a much later onset of ESD. In each section we 
mention the earlier works. Finally in section VI we summarize our findings and conclude 
with future outlook. 



2. Qubit-Qubit Interaction 

The model that we consider for our study consist of two initially entangled interacting 
qubits, labeled A and B. Each qubit can be characterized by a two-level system with an 
excited state |e) and a ground state \g). Further we assume that the qubits interacts 
independently with their respective environments. This leads to both local decoherence 
as well as loss of entanglement of the qubits. The decoherence, for instance can arise due 
to spontaneous emission from the excited states. Figure 1. show a schematic diagram 
of our model. The Hamiltonian for our model is then given by, 

n = hu;o{S% + S's) + hv{StS:^ + S+S2), (1) 

where v is the interaction between the two qubits, S^, {i =A,B) are the atomic 
energy, raising and lowering operators defined as Sf — l/2(|ei)(ei| — \gi){gi\), S^' — 
\^i){9i\ = i^rV respectively and obey angular momentum commutation algebra. We 
would use the two qubit product basis given by, 

|1) = |e)A«) |e)B \2) ^\e)A^\g)B 

|3) = b)A®|e)B \4:)^\g)A^\g)B (2) 

Now as each qubit independently interacts with its respective environment, the dynamics 
of this interaction can be treated in the general framework of master equations. The 
time evolution of the density operator p which gives us information about the dynamics 
of the system can then be evaluated from the quantum- Liouville equation of motion, 

P^~[n,p] + cp, (3) 
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Figure 1. (Color online) Schematic diagram of two qubits modelled as two two- level 
atom coupled to each other by an interaction parameter v. Here |e), \g) signifies the 
excited and ground states and luq their corresponding transition frequency. The qubits 
A and B independently interact with their respective environments (baths) which lead 
to local decoherence as well as loss in entanglement. 



where Cp includes the effect of interaction of the environment with the qubits. Note that 
in its simplest form this can be considered to be a spontaneous emission process induced 
by the vacuum fluctuation of the radiation field. For the case of simple dissipative 
environment with which the qubits are interacting independently, the effect will be 
decay of the excited state and any initial coherences of the qubit. As an example say 
for qubit A this can be written as, 

Pee = - '^lAPee 

Peg = - lAPeg- (4) 

The above equation together with the normalization Tr[p] = 1 and symmetry of the 
density matrix, define completely the dynamical system. The effect of environment as 
elucidated in equation (jl]) can be written in a compact form in terms of the atomic 
operators S*"*", S~ as, 

Cp=- ^{SlS-p-2S-pSl + pSlS-), (5) 

where the terms 7^1(75) gives the decay rate of qubit A (B) to the environment. We give 
the complete analytical solution of equation ([3]) in the basis defined by ([2]) for coupling 
to a dissipative environment ([5]) in appendix A. 

3. Concurrence Dynamics 

To investigate the effect of interaction among the two qubits on decoherence we need 
to study the dynamics of two qubit entanglement. The entanglement for any bipartite 
system is best identified by examining the concurrence [601I6T] . an entanglement measure 
that relates to the density matrix of the system p. The concurrence for two qubits is 
defined as, 

C(t) = max{0,yAr- V^- V^-^}, (6) 
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where A's are the eigenvalues of the non-hermitian matrix p(t)p{t) arranged in non- 
increasing order of magnitude. The matrix p{t) being the density matrix for the two 
qubits and the matrix p(t) is defined by, 

pW = K'^®^,f)p*W(4'^®4'^), (7) 

where p*{t) is the complex conjugation of p{t) and ay is the well known time reversal 
operator for spin half systems in quantum mechanics. Note that concurrence varies 
from C = for a separable state to C = 1 for a maximally entangled state. Though in 
general the two qubit density matrix p will have all sixteen elements, here we consider 
the initially entangled qubits to be in a mixed state [23] given by the density matrix, 

p ^l/3(a|l)(l| + d|4)(4| + (6 + c)|^)(^|); 

|^) = -=^(v^|2) + e^'^v^|3)); 
Vb + c 

^ + b + c + d ^ ^. 

where a, b, c are independent parameters governing the nature of the initial state of the 
two entangled qubits. Note that the entanglement part of the state depends on the 
initial phase x- Following ([8]) one can see that the initial two qubit density matrix have 
only six-elements. In the matrix form p is then given by. 
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(9) 



Here z = e'^^Vbc are the single photon coherences. Using the solution of the quantum- 

Liouville equation (lA.ll) it can be shown that the initial density matrix preserves its 

form for all t. Finally we calculate the concurrence defined by (El) and for the two 
qubits as, 

C{t) = Max{0,C{t)}, (10) 
where C{t) is given by, 

C{t) = 2 ||p23(t)| - VPllWP44(t)j (11) 



Let us now consider a particular class of mixed states with a single parameter a satisfying 
intially a > 0, b = c = \z\ = 1 and d = 1 — a |23]. Note that then ([8]), has the structure 
similar to a Werner state [62]. On using the dynamical evolution of the density matrix 
elements from appendix (A) and this set of initial conditions in (|TT1) . we obtain, 

C{t) = ^e-^*[(cos2x + sin2xcos2(2t;t))^/2 



a{l-a + 2w^ + w^a)], (12) 



where w = y/1 — e""^*. One can clearly see the dependence of C{t) on the interaction 
V among the qubits and the initial phase %. We see from (fT2l) that in absence of the 
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Figure 2. (Color online) Concurrence as a function of time for two initially entangled, 
interacting qubits with initial conditions 6 = c=|z| = 1.0 and two different initial 
phases x = t^/^ (black curve) and x = (I'^d curve). 

interaction v, concurrence becomes independent of the initial phase and yields the well 
established result of Yu and Eberly f23j. 
Note that C{t) can become negative if, 

a(l - a + 2w;2 + w^a) > (1 - sin^ x sm\2vt)), (13) 

in which case concurrence is zero and the qubits get disentangled. In figure (2) we show 
the time dependence of the entanglement by plotting equation ( IT2l) for f = 57 , a = 0.4 
and different values of the initial phase x- The inset of figure 2 shows the long time 
behavior of entanglement for this case. We see from figure 2 that non-interacting qubits 
(f /7 = 0) exhibit sudden death of entanglement (ESD) [visible more clearly in the inset] 
whereas when they interact (f/7 7^ 0) the concurrence oscillates between zero and non- 
zero values. Thus we see that the initially entangled qubits in presence of interaction 
V gets repeatedly disentangled and entangled leading to dark and bright periods in the 
concurrence. The magnitude of bright periods diminish with time and eventually at 
longer time this behavior vanishes completely leading to death of entanglement (ESD). 
The length of a dark period is determined by the condition (fT3|) . We have found that 
this behavior in entanglement prevails for other values of the parameter a also [59]. In 
figure 3(a) and (b) we plot the dyamical evolution of entanglement when (7 = 0), i.e 
in absence of any environmental perturbation. This is a ideal case of close quantum 
systems whose dynamics is only influenced by the initial condition of the entangled 
qubits and the inter-qubit interactions. In this case we get, 

C{t) = ^ [v'cos2x + sin2xcos(2t;t) - V^a(l - a)] (14) 

For both the case of a = 0.2 and a = 0.4 with an initial phase of x/4 , we observe 
sinusoidal behavior of entanglement as seen in both (a) and (b) of figure 4. Thus 
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Figure 3. (Color online) Evolution of Concurrence for two initially entangled, 
interacting qubits with initial conditions a — 0.4, b = c — \z\ — 1.0 and different initial 
phases x- Here 7 = 0. The magnitude of bright periods in absence of environment 
does not diminish in magnitude. 



there is no ESD in absence of the environment in this case. For another value of the 
initial phase x/2 we observe dark and bright periods of entanglement. The periods of 
disentanglement (dark periods) are governed by the condition ^Ja{l — a) > \ cos(2t>t)|. 
It is clearly visible from the plots that in absence of any environment the amplitude of 
the bright periods does not diminish at all and thus the qubits gets back their initial 
entanglement completely. This regeneration of entanglement is due to the inter-qubit 
interactions. Note that similar behavior in concurrence dynamics (collapse and revival 
of entanglement) have been predicted in earlier studies of non-interacting qubits in 
atom cavity systems. For example it was shown that for double Jaynes-Cumming (JC) 
[63] model, with completely undamped non-interacting cavities entanglement shows a 
periodic death and re- birth feature [M]- This was attributed to exchange of information 
between the finite number of cavity modes and the atoms - a new kind of temporary 
decoherence mechanism. In another work pairwise concurrence was calculated among 
four qubits, where the qubits were formed by the cavity modes and atoms [65]. Here 
again JC like interaction between the atom and cavity gives rise to dark and bright 
period in the entanglement dynamics of the qubits. It was shown that during the 
period when the concurrence between the cavities vanish, the concurrence between the 
atoms reaches its peak and vice- versa. This only happens as the cavities where assumed 
to be lossless with finite number of mode and thus without environmental decoherence. 
Further it was shown that for qubits remotely located and in contact with their respective 
environment when driven independently by single mode quantized field, one gets dark 
and bright periods of entanglement instead of ESD, a feature similar to single atom 
behavior in cavity quantum electrodynamics [66], [67] . These works [63], [611 [HSl [66], [67] 
differ from our's as we focus on the effect of interaction among the qubits in presence of 
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a decohering environment. Note that in a more recent work it was shown how oscillators 
interacting with a correlated finite temperature Markovian bath can lead to dark and 
bright periods in entanglement for certain initial conditions [19]. 

4. Pure Dephasing of the Qubits 




Figure 4. (Color online) Schematic diagram of two qubits modelled as two two-level 
atom coupled to each other by an interaction parameter v. Here |e), \g) signifies the 
excited and ground states and loq their corresponding transition frequency. The qubits 
A and B independently dephase to their respective environments (baths) which leads 
to decoherence and thus loss in entanglement. The corresponding dephasing rates are 
given by Ta and Tb respectively. 

In order to demonstrate the generic nature of our results, we consider other 
models of the environment. A model which has been successfully used in experiments 
[68] involves pure dephasing. The mathematical formulation for this kind of an 
environmental model can be done via a master equation technique and is given by, 

Cp=- US^S^P - '^StpSt + pStSt) (15) 

i=A,B 

where Va{Xb) is the dephasing rate of qubit A (B). Substituting ( fTSj) in ([3]) we get 
the equation for dynamical evolution of the qubits under the infiuence of this kind of 
an environment. Note that in this model the populations do not decay as a result 
of the interaction with the environment whereas the coherences like p2z{t) decay as 
P23(0)e~*^^'*"'"'"-^^*. Let us now study the the effect of interaction v between the qubits on 
the dynamics of entanglement. We assume the same initial density matrix of equation 
([9]) with the initial conditions d = 1 — a,6 = c= \z\ = 1 and a > to calculate the 
concurrence. One can clearly see from the solution of quantum-Louiville equation given 
in appendix (B) that under pure dephasing, the form of matrix in ([9]) is preserved for 
all time. Using the solutions of the master equation ([3]) derived in appendix (B) for the 
environment effects given by ( ITSl) and substituting in equations ( |TOl) . ( ITTi) we get the 
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Figure 5. (Color online) Concurrence as a function of time with initial conditions 
b — c — \z\ = 1 and two different values of the initial phase x ^^r the dephasing 
model. The red and black curve in figure is for % = 7r/4 and tt/2 respectively. Here 
the interaction parameter is taken to he v/T = 4. 



time dependent concurrence for this model to be, 

Coit) = '^[e-^{e-^^cos^x + sin\{cos{QiT) 

-l-sm{n,T)ry/'-^a{l~a)], (16) 

where the suffix D signifies that the concurrence is calculated for a dephasing 
environment and we assume F^i = = F. Here t = Tt and fli = ^ (2f /F)^ — 1. 

For f = we get Coit) = 2/3 [e^^"^ — — a)], which is independent of the initial 

phase X- We find death of entanglement for r > (1/2) ln[l/y/a(l — a)]. Note that Yu 
and Eberly [25] have considered this case earlier but for a = 1 only, in which case there 
is no ESD. In figure (5) we show the time dependence of entanglement for a purely 
dephasing model, for a = 0.2 and initial coherences governed by the phase x- From the 
figure we see that for v 0, the two qubit entanglement exhibits the phenomenon of 
dark and bright periods. Further we also see that for v 0, dark and bright periods 
continues beyond the time when ESD occurs for noninteracting qubits. This kind of 
behavior in the entanglement dynamics is found for other values of the parameter a and 
interested readers are refered to [59] for further dicsussions on these. 
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5. Concurrence Dynamics in Correlated Environmental Models 
5.1. Effect of Correlated Dissipative Environment 

We next consider an environment involving correlated decay and show how coupling to 
such environment can lead to new effects in the entanglement dynamics for two qubit 
systems. We will consider the case of both non-interacting as well as interacting qubits 
for this model of the environment. To keep the analysis simple and get a better physical 
insight on the decoherence effect of this environment we will first study the case of 
non-interacting qubits. We assume as before that the qubits interacts independently 
with their respective environments with decay rates 7a and 75. Further we assume 
that the qubits are close enough (r << A, r being the inter-qubit distance and A the 
wavelength of emitted radiation in process of a decay) such that they can undergo a 
correlated decay with decay rates Tab(Xba) for qubits A{B). Whether this would lead 
to further decoherence is a question we want to investigate. Note that the entanglement 
dynamics of two non-interacting two level atoms in presence of dissipation caused by 
spontaneous emission was studied earlier in details by Jakobczyk and Jamroz [37] . 
They even considered correlated model of dissipative environment and showed possible 
destruction of initial entanglement and possible creation of a transient entanglement 
between the atoms. Further they also discussed the question of non-locality and how 
it is influenced by the spontaneous emission by explicitly showing the violation of Bell- 
CSHS inequality. One of the chief difference between this work and ours is the initial 
density matrix p considered and the interaction introduced between the qubits . While 
we consider the possibility of both the qubits (atoms) to be initially excited and show 
its important consequences on the decay dynamics, they have neglected this effect by 
putting pii(O) = 0. We would show later in this paper (as can also be seen from 
their results) that the dissipative environment preserve the form of the initial p. Hence 
Pii{t) = for all time in their case. Moreover, in a recent work the entanglement 
dynamics of two initially entangled qubits for collective decay model was studied in 
context to ESD, by Ficek and Tanas [38]. They considered an initial density matrix of 
the from. 



It can be clearly seen that in this case the two-qubits are initially prepared in an 
entangled state by the two-photon coherences. They further show that for this initial 
condition the single photon coherences are never generated. Moreover the dipole-dipole 
interaction that they consider for the two qubit system have no influence for this 
initial condition. Ficek and Tanas predicted dark periods and revival in the two qubit 
entanglement in their work due to the correlated nature of the bath, we on the other hand 
consider the initial density matrix of the form with single photon coherences and 
show that any coherent interaction among the qubits does influence the entanglement 
dynamics at all later time. 



P 




(17) 
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We now include the effect of a dissipative environment with both independent and 
correlated decay of the qubits via a master equation technique given by, 

j,k=A,B ^ 

= 7, (18) 

The time evolution of the density operator p which gives us information about the 
dynamics of the system can then be evaluated by solving the quantum-Liouville equation 
([3]) with the environmental effect included by equation f|T8|) and taking v = 0. Next as 
before we consider the qubits to be intially entangled with their initial state to be a mixed 
state defined by the density matrix ([9]). We then solve the quantum-Louiville equation 
to study the dynamical evolution of the system. The reader is refered to appendix (C) 
for explicit solution of the time dependent density matrix elements. One can clearly see 
from appendix (C) that for this kind of model of the environment, as before the initial 
density matrix preserves its form for all time t. Now using appendix (C) in equations 
( ITOj) and (ITT]) and the initial conditions a > 0, d = 1 — a, b = c = \z\ = 1, we obtain the 
concurrence dynamics of two initially entangled non-interacting qubits for this model of 
the environment as, 

C{t) = -e-''*{[{cosxcosh(n) - sinh(n) + aC(t)}^ 
3 



+ sin2x]^/'- V3a[l-«:(t)]}, 

(19) 



where ((t) and K(t) are given by 




(20) 




^ (tTTT^) " ^ V-*{cosh(n) 

- cosxsinh(rt)}, (21) 

For simplicity we have assumed equal decay rates of both the qubits, 7a = 7_b = 7 
and Tab = Tba = T- One can clearly see the dependence of C(t) on the correlated 
environmental effect given by F and the initial phase x equation f|T9|) . We see from 
( IT9|) . ( l20l) and ( 12T1) that for F = 0, concurrence becomes independent of the initial phase 
and yields the result of Yu and Eberly [23j. Note that C{t) can become negative if, 

3a[l - K{t)] > [{cosxcosh(Ft) - sinh(Ft) + aC(t)}^ 

+ sin^ x] (22) 
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Figure 6. (Color online) Time evolution of concurrence for a = 0.2, b = c = \z\ = 1 
and two different initial phases x for two non-interacting qubits in contact with a 
correlated dissipative environment. Here r/7 = signifies absence of common bath 
for the qubits. 
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Figure 7. (Color online) Time evolution of concurrence govern by the initial condition 
a — 0.4 for two non-interacting qubits in contact with a correlated dissipative 
environment. Here all other initial parameters remains the same as fig (7). r/7 = 
signifies absence of any common bath in which case entanglement sudden death (ESD) 
is observed. 



in which case concurrence is zero and the qubits get disentangled. To understand how 
correlated decay of the qubits might effect their entanglement we study the analytical 
result of equation (fT9|) for different values of the parameter a and x- In figure (6) we show 
the time dependence of entanglement for a = 0.2 and two different values of initial phase 
X and correlated decay rate of F = O.87. Note that for F = 0, there is no ESD in this 
case [23] and concurrence monotonically goes to zero as t — > 00. For F 7^ we observe 
new behavior in the entanglement of the qubits. Concurrence is seen to have a much 
slower decay in comparison to when F = 0. For a initial phase of x = 7r/4 we observe 
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that the condition in equation fl22p is satisfied and entanglement vanishes temporarily 
i.e the qubits get disentangled. The entanglement gets regenerated at some later time 
and finally goes to zero very slowly as t — > oo. Note that this disentanglement and 
re-entanglement phenomenon is non periodic and is very sensitive to initial coherence 
among the qubits, for example it do not occur when the initial coherence is governed 
by the phase x = ^/2- In figure (7) we plot concurrence for a = 0.4. For this value of a 
ESD is observed for F = but not for F 7^ 0. Instead we observe disentanglement and 
regeneration of entanglement among the qubits for % = 7r/4. Here again we find that no 
dark and bright periods nor any ESD for initial phase of x = 7r/2. Further, note that 
for initial phase x = we have a longer time interval during which the qubits remain 
disentangled before getting entangled again, in comparison to the case for a = 0.2. Thus 
we find that the time interval between disentanglement and regeneration of entanglement 
as well as the magnitude of regeneration strongly depends on the initial coherences of 
the initially entangled qubits. Hence we can conclude that for non-interacting qubits in 
contact with a dissipative correlated environment no ESD occurs. 

Let us now consider the case of two initially entangled interacting qubits in contact 
with the correlated environment. The dynamical evolution of the system in presence 
of interaction v for correlated model of environment is evaluate in details in appendix 
(D). We use the solutions of appendix (D) in (fTT!) to calculate the concurrence for this 
environment. Note that the solutions are essentially valid under the assumption that 
our initial two qubit density matrix p is given by equation ([9]). Further we consider 
as before that the two entangled qubit's evolution is governed by the initial conditions 
a > 0,d = 1 — a, b = c = \z\ = 1. Hence the time dependent concurrence for two 
initially entangled interacting qubits becomes, 

C{t) = ^e"^*{[{cosxcosh(Ft) - sinh(Ft) + aC(t)}^ 

+ cos2(2ft) sin^ xY^^ - v^3a[l - ^(t)]}, 

(23) 

C{t) = Max{0,C{t)} ; (24) 

where ({t) and K,{t) are given by equations fl20l) and fl2Tl) respectively. The dependence 
of concurrence for C{t) > on the interaction strength v between the qubits is clearly 
visible in equation (1231) . Further now we can see that the condition of complete 
disentanglement of the qubits is given by, 

3a[l - K{t)] > [{cosxcosh(Ft) - sinh(Ft) + aC(t)}^ 

+ cos2(2ft)sin2x] (25) 

When condition fl25l) is satisfied, C(t) < and hence C{t) = 0. Next to study the effect 
of qubit-qubit interaction on the entanglement dynamics we plot the time dependent 
concurrence for different value of a, initial phase x and correlated decay rates of F = O.87 
in figures (8) and (9). We observe in the figures that for an initial phase of x = 7r/2 
concurrence exhibits dark and bright periods at initial time for both a = 0.2 and 
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Figure 8. (Color online) Time evolution of concurrence for interacting qubits 
in contact with a correlated dissipative environment with correlated decay rate of 
r/7 = 0.8. Here b = c = \z\ = 1. A long period of disentanglement is observed for 
initial phase x = 7r/4. Here the interaction strength among the qubits is taken to be 
w/7 — 5.0 




Figure 9. (Color online) Time evolution of concurrence for interacting qubits in 
contact with a correlated dissipative environment for same parameters as figure (9) 
and a — 0.4. The dark and bright periodic features sustain for a longer time for initial 
phase of X = 7r/2. Much longer period of disentanglement now observed for x — 7r/4. 



a = 0.4. For longer time the concurrence shows a damped oscillatory behavior. We 
attribute this effect to the competition between the fast inter-qubit interactions v and 
the environmental decays. For longer time the correlated decay becomes dominant 
and leads to a slow damped oscillatory decay of the entanglement. For x = 7r/4 the 
dark and bright periods are not very pronounced and is over shadowed very quickly by 
the correlated decay. Note that for this value of initial phase we find that there exist 
a long period of time during which the qubits remain disentangled. At a much later 
time entanglement gets regenerated and increases initially and then starts decaying very 
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slowly after . This behavior is quite different from the dark and bright periods seen for 
other models of the environment. Thus we see that for interacting qubits there is no 
ESD for this model of the environment. Instead we find dark and bright periods with 
long period of disentanglement whose occurrence depends on the initial coherence. 

5.2. Delay of ESD by Correlated Dephasing Environment 




Figure 10. (Color online) Schematic diagram of two qubits modelled as two two- 
level atoms. Here Wg is the transition frequency of the excited state |e) to the ground 
state \g). The qubits A and B independently dephase to their environments (baths) 
with a dephasing rate of Ta,^b respectively. The qubits can also interact with the 
environment collectively when they are at proximity giving rise to correlated dephasing 
represented by the decay rate Fq. 



Finally we consider a purely correlated dephasing model of the environment and 
study the effect of such an environment on the entanglement dynamics of two qubits. 
Note that this kind of model is popular among solid state systems like semiconductor 
quantum dots. We will study the behavior of entanglement for both non-interacting and 
interacting qubits. As before to keep our analysis simple and to get a better physical 
insight to the question of decoherence for this kind of environment we will first study 
the case of non-interacting qubits. We will then generalize our results by introducing 
the interaction among the qubits. For non-interacting qubits the Hamiltonian for our 
model is given by ([T]) with f = 0. The effect of the dephasing environment on the qubits 
is included via a master equation technique and is given by, 

Cp= - T,{StStp-2StpSt + pStSt) 

i=A,B 

— 2Tq{S\S%p — S^pS^ + pS\S% — S^pS^), 

(26) 

where TAi^s) and 2ro are respectively the independent and correlated dephasing rate 
of qubit A (B) . The dynamical evolution of this system can then be studied by solving 
the quantum-Louiville equation ([3]) for v = and including the effect of environment by 
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using (I26p . We now consider as earlier that the initial state of the two qubits is defined 
by the density matrix p ([9]). Then the solution of the quantum-Louiville equation for 
this model of the environment is given by, 

\c, (27) 

(28) 

P22(t)-P33(t) (29) 

All other matrix elements of the two qubit density matrix p are zero. Now using the 
solutions of (l29l) it is straight forward to show that, for pure dephasing of the qubits, the 
form of matrix in ([9]) is preserved for all time. Note that in such a model the populations 
do not decay as a result of the interaction with the environment whereas the coherences 
like P2'i{t) decay as ~ p2z{^)e~^^^^^^~'^^°'>^ for x = or mod vr. Let us now study the 
effect of correlated dephasing of the qubits on the dynamics of entanglement. For the 
initial conditions d = 1 — a^h = c = \z\ = 1 and a > on using (l27|l in (fTTj) we get the 
expression for time dependent concurrence as, 

Cz.(t) = ^{e-^(^-^°)*-y^(r^} (30) 

where we have assumed = = F for simplicity. From equation (l30l) it is clearly seen 



Pll{t) = gfl, P22{t) = -b, P33(t) = 

P23(t) = ^k|e-(^-+^--^^°)V>^ 

P32{t) = Plsit), pii{t) = 1 - pii(t) - 




Figure 11. (Color online) Time evolution of concurrence for two non-interacting 
qubits in contact with a purely dephasing environment for initial condition given by 
a = 0.2 and 6 = c = \z\ = 1. The effect of correlated dephasing shows up as delay in 
the onset of ESD. 



that in a purely dephasing environment entanglement among the qubits is independent 
of the initial coherence given by x aiid depends only on a and Fq. In figure (11) we plot 
the time dependence of concurrence for a = 0.2. We find that the effect of correlated 
dephasing is manifested in the delay of the onset of ESD. The time for the onset of 
ESD is given by t > 1/2(F — Fo){l/ln -y/a(l — a)}. From the figure its is clealy visible 
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that with increase in correlated decay Fq, the onset of ESD gets delayed further until 
To = r, when concurrence bceomes independent of the dephasing rates and is given by 
C = 2/3[l — ^Ja{l — a)]. This situation represents a decoherence free subspace where 
concurrence becomes solely dependent on the value of a i.e population of the excited 
state of the two qubits. Note that this kind of situation has already been tailored to 
study entanglement in decoherence free subspace [68j. 

Let us now include the interaction among the qubits and study how this interaction 
might influence the entanglement dynamics for this model of the environment. The 
Hamiltonian of the two qubit system and its coupling to the environment is then given 
by equations and (|26ll respectively. To study the dynamics of entanglement we follow 
a similar process as described earlier. We use the solution of quantum-Louiville equation 
derived explicitly in appendix (E) and substitute them in equation (fTTl) to calculate the 
time dependence of concurrence C. With the initial conditions a = 1—d, b = c = \z\ = 1, 
then we get, 

Co{t) = l[e-^^-'^^'{e-'^^~^^^'cos'x 

+ sm\{cosin't)-^^^^^smin't)fY/^ 

- ^a{l-a)]; (31) 
C(t) = Max{0,C'fl(t)} (32) 

where fl' = ^iv"^ — (F — Fq) and we have assumed F^ = F^. One can clearly see 
the dependence of concurrence on the interaction v among the qubits for Cd > 0. Note 
that due to the interaction among the qubits now concurrence becomes dependent of the 
initial phase x- To understand the behavior of entanglement in presence of interaction 
{v/'y = 5.0) among the qubits we plot the time dependence of concurrence for different 
initial phase x ci-nd correlated dephasing rates Fq in figures (12-14). We consider the 
case, a = 0.2 only to do a comparative study on the behavior of concurrence in presence 
and absence of inter-qubit interactions. Note that we have already discussed the effect 
of correlated dephasing on the two qubit entanglement for this value of a. Let us now 
focus on any new feature that arises due to the qubit-qubit interactions. We can see 
clearly from figure (12) that for f ^ 0,x = 7r/4 the two qubit concurrence shows a 
damped oscillatory behavior which leads to dark and bright periods at longer time 
before eventual death of entanglement. The generation of dark and bright periods is 
seen to delay the death of entanglement even further in comparison to that induced by 
correlated dephasing in absence of qubit-qubit interactions. Moreover in figure (13) we 
see that both the oscillatory behavior as well as dark and bright periods is enhanced 
with an increase in correlated dephasing rate. When we change the initial phase to 
7r/2 for Fq = 0.2F we find (figure 14) no oscillatory behavior in entanglement rather 
a completely dark and bright periodic feature with eventual delayed death. Thus we 
see that the onset of dark and bright periods for this kind of environment model is 
profoundly influenced by the initial coherence of the two qubit system. 
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Figure 12. (Color onlinc)Timc evolution of concurrence for two interacting qubits 
with interaction strength v/T = 5.0 in contact with a purely correlated dephasing 
environment and initial condition a — 0.2. ,x = 7^/4:, b = c = \z\ = 1. The red 
curve correspond to concurrence of non-interacting qubits. Concurrence is seen to 
exhibit initial oscillations followed by dark and bright periods with eventual death of 
entanglement in presence of interaction. The interaction also leads to delayed death 
of entanglement. Here Tq is the correlated dephasing rate. 
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Figure 13. (Color online)Time evolution of concurrence for two interacting qubits in 
contact with a correlated dephasing environment with same parameters as for figure 
(12) but higher correlated dephasing rates. The effect of higher correlated dephasing 
manifests itself by increasing the periodicity of dark and bright features in concurrence 
. Here again we find that dark and bright periods is followed by death of entanglement. 

The phenomenon of dark and bright periods in entanglement should have direct 
consequences for systems like ion traps , quantum dots, the later being currently the 
forerunner in implementation of quantum logic gates. The interaction between qubits 
considered in this paper are inherently present in these systems. In quantum dots for 
example, ~ few ns and one can get a very large range of the parameter F"^ (1-100's 
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Figure 14. (Color online)Tiine evolution of concurrence for two interacting qubits in 
contact with a correlated dephasing environment with initial condition a = 0.2. ,x = 
7r/2, 6 = c — \z\ — 1. Concurrence is seen to be sensitive to initial coherence among 
the two qubits. It does not exhibit initial oscillations for this value of x but dark and 
bright periods with eventual death of entanglement in presence of interaction. Here 
the interaction strength is taken to be v/T = 5.0 

of ps) [69] • Further the interaction strength v can have a range between l/xev - 1 mev 
depending on gate biasing [131 [TOl [7T1 [72] . An earher study [73] reports 7 ~ 40 — lOO/xev 
and couphng strength of ~ 100 — 400yuev, thereby making f/7~l — lOfor quantum dot 
molecules. Thus experimental parameters are in the range we used for our numerical 
calculation. 

6. Summary 

In summary we have done a detail study of decoherence effect for non-interacting 
and interacting initially entangled qubits in contact with different environments at 
zero temperature . We have shown how the interaction between qubits generates the 
phenomenon of dark and bright periods in the entanglement dynamics of an intially 
entangled two qubit system in contact with different environments. We found this 
feature of dark and bright periods to be generic and occurs for various models of 
the environment , as an example in a correlated dissipative environment we found 
the phenomenon of dark and bright periods in entanglement dynamics even though 
there is no sudden death of entanglement. Moreover for purely dephasing models of 
the environment we found that the dark and bright periods feature sustains longer and 
delays the sudden death. We found that there is no sudden death of entanglement for a 
correlated dissipative environment but rather depending on the initial coherences in the 
system entanglement can show a substantial slower decay and even the phenomenon 
of dark and bright periods. For a simple pure dephasing environment as well as for 
correlated dephasing environment we have shown the existence of sudden death of 
entanglement. Due to correlated dephasing we found delayed death of entanglement. 
Further, in the correlated dephasing model we found that the onset of dark and bright 
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periods is sensitive to the initial coherence in the system. The frequency of dark and 
bright periods was found to depend on the strength of interaction between the qubits 
as well as on the correlated decay and dephasing rates. As a future perspective it would 
be interesting to study the effect of qubit-qubit interaction for environments having 
temperature fluctuations. Further it would also be interesting to extend our study to 
multi-qubit entanglement. An important class of states that can be treated for this 
purpose are the GHZ and W states. Moreover cluster states [74J can also be considered 
as other probable candidates for the study of decoherence and loss of entanglement. 
This work was supported by NSF grant no CCF-0829860. 



Appendix A. Solution of the quantum-Louiville equation in the two qubit 
product basis for two interacting qubits in contact with a dissipative 
environment. 



P22{t) 

P33W 

Pl2(t) 
Pl3(t) 



^P22(0)e-^*(1 + cos{2vt)) + lp33(0)e-^*(l - cos(2^;t)) 



+ Pii(0)e 



l)-o(P32(0)-p23(0))e-^*sin(2t;t) 



+ 



+ 



P24(t) 



^P22(0)e-^*(1 - cos(2^;t)) + ip33(0)e-^*(l + cos(2t;t)) 

pii(0)e-2^*(e^* - 1) + ^(P32(0) - P23(0))e-^* sm{2vt) 

Pi2(0)e-3T*/2 cos{vt) + 2pi3(0)e-3T*/2 sm{vt) 
Pi3(0)e-3^*/2 pog(^^) ^ zpi2(0)e-3^*/2 gj^^^^) 

Pi4(0)e-^* 

^(P22(0) - P33(0)) sm{2vt)e'^' + ^p23(0)e-^*(l + cos{2vt)) 
^P32(0)e-^*(l-cos(2t;t)) 



P24(0)e-^*/2 cos{vt) - ^p34(0)e-^*/2 ^-^^^^^ _ ^^^^q^ f , , 



y + (9/4)72 



P34(i) 



+ e-^*/^{2wcos(wt) - 3i^sm{vt)}] - pi3(0) 

- e~^*/2{2wsin(wt) + 37cos(t;t)}] 

= P34(0)e-^*/2 cos{vt) - ip24(0)e-^*/^ sin(t;t) 

- Pl2(0) 



+ (9/4)72 



[376 



-2-ft 



1 



- Pl3(0) 



[876-2^* - e~^^/^{2vsm{vt) + 37cos(wt)}] 

[2it;e-2^* + e-^^/^{2v cos{vt) - 8^7 sin(t;t)}] 

(A.l) 
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and p32(t) = phit),P2i{t) = pl^{t),p3i{t) = pl^{t),pii{t) = Pl^{t),pi2{t) = 
P24(^)'P43(^) = Plii^), Pu{t) = 1 - Pii(t) - P22(i) - p33{t)- Note that here we have 
considered Ja = Ib = 1- 



Appendix B. Solution of the quantum-Louiville equation in the two qubit 
product basis for two interacting qubits in contact with a purely dephasing 
environment. The solutions correspond to the initial matrix p defined in 
equation ([9]) 



Pii(t)=Pii(0) 

P22{t) = ^P22(0) 



+ 2^33(0) 



1 + e-(r^+r^)*/2 f (2fit) + ^^^^ sin {2nt) ] 



{B.l] 



An 



1 _ ^HrA+rsW2 I^Qg (2fit) + ^^^^J""^ sin (2fit) 



^.[P23(0)-P32(0)K-(r.^r.)V2^^^^^^^^ 



2n 



P33{t) = -P22(0) 
+ ^P33(0) 



1 _ e-^r,+r,)t/2 f (2fit) + i£A±££ 



(B.2) 



An 



■ sin {2nt) 



1 + e-(r^+rB)i/2 L^g ^2nt) + ^^^^"^^^^ sin (2^^) 



«[P23(0) - p32i0)]ve 



'irA+rB)t/2 



2n 



P23(t) = ie-(^-+^-)*/2[p23(0){e 



■ sin {2nt) 



<r,+vM2 ^ (2fit) + ^^^^-^ si 



An 

+ P32(0){e"(r-+r-)*/2 _ (2f^t) - ^^^^^^ sin (2^]^)}] 
+ sin {2nt) [p22(0) - P33(0)], 



(B.3) 
sin {2nt)] 



n 



P32 (t) = p*23 (t) , P44 (t) = 1- pii (t) - P22 (t) - P33 (t) . 



(B.4) 
(B.5) 
(B.6) 



All other elements of the density matrix p defined in the two qubit product basis 
remains zero for all time t. 
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Appendix C. Solution of the quantum-Louiville equation in the two qubit 
product basis for two non-interacting qubits in contact with a dissipative 
environment which results in correlated decay. The solutions correspond to 
the initial matrix p defined in equation ([9]) 



Pii(t)=Pii(0)e-2^* ; 

P22(t) = ^P22(0)e-^*(1 + COsh(ri2t)) - ^P33(0)e-^*(1 - COsh(ri2t)) 



(C.l) 



- Pii(0)e 



-2jt 



7' + r 



1 



' 1 + ri2^ 



[p23(0)+P32(0)]e-^*sinh(ri2t) 



.7-r 



^7-r 



12 



-Fiat 



12 , 



.7 + r 



12 , 



(C.2) 



P33(i) = ^P33(0)e-^*(1 + COsh(ri2t)) - ^P22(0)e-^*(1 - COsh(ri2t)) 



- Pii(0)e 



-27* 



7' + r 



12 



^2 _ r2 

,7 -Ll2, 



^[p23(0)+P32(0)]e-^*sinh(ri2t) 



'7 + r 



v7 



12 \ ^r,,* ^ 



12 . 



7 



12 



.7 + r 



12 . 



(C.3) 



P23(t) = ^P23(0)e-^*(1 + COsh(ri2t)) - ^P32(0)e-^*(1 - COsh(ri2t)) 



:[P22(0) +P33(0)]e-^* sinh(ri2t) - Pii(0)e- 



-2-it 



12 



2 _ p2 



.7 



12, 



+ 2^11(0)^"'* 



12 



v7 - ri2. 



'7-r 



12 



.7 + ri2. 



P32(t)=P23W> P44(t) 



Pll(t) -P22(t) -P33(t)- 



(C.4) 
(C.5) 



All other elements of the density matrix p defined in the two qubit product basis 
remains zero for all time t. 



Appendix D. Solution of the quantum-Louiville equation in the two qubit 
product basis for two interacting qubits in contact with a dissipative 
environment which results in correlated decay. The solutions correspond to 
the initial matrix p defined in equation (I9l) 



Pii(t)=Pii(0)e-2^* 



(D.l) 



P22(t) = ip22(0)e-^*(cos(2t;t) + cosh(ri2t)) - ^p33(0)e-^*(cos(2t;t) - cosh(ri2t)) 



-pii(0)e-2^*'^ 



.7 



12 , 



12 



7 



12 



-ri2t 



.7 



12 . 



.7 + ri2. 
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- ^[P23(0) +P32(0)]e-^*sinh(ri2t) + ^[p23(0) -P32(0)]e-^*sin(2t;t); 



(D.2) 



^Ps3{0)e'^\cos{2vt) + cosh(ri2t)) - ^P22(0)e^*(cos(2i;t) - cosh(ri2t)) 



Pii(0)e-2^* 
1 



7' + r^ 



12 



+ oPii(0)e 



-7t 



12 , 



12 



.7-r 



J^izt _|_ I J 1 ^-il2 



12 . 



7 + r 



12 \ ^_ri2t 



12 . 



2[P23(0) +P32(0)]e^^*sinh(ri2t) - -[p23(0) - P32(0)]e-^* sin(2t;t); 



(D.3) 



P23W = ^P23(0)e-^*(cos(2t;t) + cosh(ri2t)) - ip32(0)e-^*(cos(2t;t) - cosh(ri2t)) 



Pii(0)e 



-2'yt 



12 



12, 



V7' - 

1. 1 



't + Ti 

.7-r 



+ 



^7-r 



12 



12 , 



7 + r 



12 , 



-[P22(0) +P33(0)]e-^*sinh(ri2t) + -e-^*[p22(0) -P33(0)]2sin(2^;t); 



(D.4) 



P32 (t) = p^g (t) , P44 (t) = 1 - Pll (t) - P22 (t) - P33 (t) . 



(D.5) 



All other elements of the density matrix p defined in the two qubit product basis ([3]) 
remains zero for all time t. 

Appendix E. Solution of the quantum-Louiville equation in the two qubit 
product basis for two interacting qubits in contact with a purely dephasing 
environment. The solutions correspond to the initial matrix p defined in 
equation (I9D 



Pii(t)=Pii(0) 

P22(t) = ^P22(0) 



2P33(0) 



(E.l) 



1 + e-(rx+rB-2ro)t/2 ^2fi't) + ^ sin (2fi't) 

1 - e-(r-4+r,-2ro)V2 (2fi't) + + ~ sin (21]'t) 



+ -[P^3(0)-P32(0)K-(r..^r.-2ro)V2^^^^^^^^^ 



2^]' 



(E.2) 



P33W = 2^22(0) 
1 



1 - e-(r-4+rB-2ro)V2 (2fi't) + ^ sin (21]'t) 

1 + g-(r,+r,-2ro)t/2 ^2fi't) + + ~ sin (2n't) 



+ 2^33(0) 

_ ^[p23(0)-p32(0)]t;e-(^^+^^-^roW2 

2^]' 



sin (2^]'^) 



(E.3) 



P23(t) = ie-(r-+^--2r«)*/2[p23(0){e-(^-+r--2^°)*/2 
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+ cos i2n't) + + ~ sin (2fi't)} 

+ P32(0){e-(r-+r--2r")*/2 _ (2^]'^) - + ~ sin (2r]'t)}] 

^^g-(rA+rB-2ro)t/2 
+ sin (2^]'^) [p22(0) - P33(0)] , 

(E 



P32(t) = P*2s{t), pu{t) = 1 - pn(t) - P22(t) - P33(t). (E.6) 

All other elements of the density matrix p defined in the two qubit product basis ([3]) 
remains zero for all time t. 
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